Exact-exchange kernel of time-dependent density functional theory: Frequency 
dependence and photoabsorption spectra of atoms 
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In this work we have calculated excitation energies and photoionization cross sections of Be 
and Ne in the exact-exchange (EXX) approximation of time- dependent density functional theory 
(TDDFT) . The main focus has been on the frequency dependence of the EXX kernel and on how it 
affects the spectrum as compared to the corresponding adiabatic approximation. We show that for 
some discrete excitation energies the frequency dependence is essential to reproduce the results of 
time-dependent Hartree-Fock theory. Unfortunately, we have found that the EXX approximation 
breaks down completely at higher energies, producing a response function with the wrong analytic 
structure and making inner-shell excitations disappear from the calculated spectra. We have traced 
this failure to the existence of vanishing eigenvalues of the Kohn-Sham non-interacting response 
function. Based on the adiabatic TDDFT formalism we propose a new way of deriving the Fano 
parameters of autoionizing resonances. 

PACS numbers: 31.15.Ew, 31.25.-v, 71.15.-m 
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I. INTRODUCTION 

The photoabsorption spectrum of an atom or molecule 
reveals detailed information about the structure and dy- 
namics of its constituent electrons. Important phenom- 
ena, such as Fano resonances^ and multiple-particle ex- 
citations, which can solely be attributed to the electron- 
electron interaction, has been under extensive experimen- 
tal and theoretical investigation for several decades. To 
describe and quantify these exciting many-body effects 
is, however, still a major theoretical challenge. Tradi- 
tional approaches to this problem are based on infinite- 
order many-body perturbation expansions which lead to 
complicated integral equations already at the Hartree- 
Fock level of approximation. This level of approxima- 
tion is also referred to as the random phase approxima- 
tion with exchange (RPAE) and was extensively used by 
Amusia^'^ and by Wendin"* for the description of atomic 
photoabsorption spectra. From a mathematical or nu- 
merical point of view this entails the solution of an inte- 
gral equation known as the Bethe-Salpeter equation for a 
four-point vertex function, the solution of which is com- 
putationally very demanding in low-symmetry systems 
like larger molecules or nano structures. 

Since the pioneering work of Andc^ and Zangvill and 
Soven^ in the late seventies, time-dependent density func- 
tional theorj«li^ (TDDFT) has emerged as a more con- 
venient route toward the calculation of spectral prop- 
erties. Within TDDFT, as it is presently formulated, 
only two-point correlation functions appear and, in ad- 
dition, one arrives at explicit expressions for the lin- 
ear density response function, thus avoiding the need 
for integral equations. The price to pay is the diffi- 
culties associated with obtaining approximations to the 
frequency-dependent exchange-correlation (XC) kernel, 
constituting one of the basic elements of TDDFT. In 
most calculations until today one uses a static, frequency- 
independent, kernel. These so called adiabatic approxi- 



mations have, in many cases, produced relatively accu- 
rate and important results - especially for static prop- 
erties or small excitation energies. There are, however, 
many cases where the inadequacies of the adiabatic ap- 
proximations are easily discernible. 

Until recently, no systematic way of obtaining succes- 
sively better approximations to the XC kernel of TDDFT 
has been available. By means of the variational formu- 
lation of many-body theorj*^ we have recently developed 
a tool for constructing approximations to the XC ker- 
nel, which not only incorporates the important energy 
dependence of the kernel but does it in a way so as to 
guarantee the fulfillment of many sum rules and conser- 
vation lawsji^ As discussed in these previous publications 
the variational approach is not uniquely defined and some 
implementations are preferable to others. For instance, 
starting from the functional of Luttinger and Wardii 
(LW) is usually a better approach as compared to start- 
ing from the Klein functional. The simplest approxi- 
mate kernel derived from the latter functional is that of 
the time-dependent exact-exchange (TDEXX) approxi- 
mation. This approximation can be considered as an 
attempt to mimic the time-dependent Hartree-Fock ap- 
proximation within a TDDFT framework, a point which 
has been discussed at length in previous workJ^ii^ The 
exact-exchange (EXX) kernel has, however, yet received 
little attention especially with regard to its ability to de- 
scribe spectral properties which necessitates a calcula- 
tion on the real frequency axis. Previous investigations 
have either addressed the adiabatic limit, ^^'^^ focused on 
extended systems at low frequenciesii or most recently 
carried out model calculations in the time domainJS 

The present paper presents a numerical and analytical 
investigation of the fully frequency-dependent EXX ker- 
nel of closed-shell atoms. It is the fourth and latest in 
a series of papers by us concerned with approximations 
obtained from the variational scheme to TDDFT. In the 
first paper— we calculated the EXX kernel of the electron 
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gas. From that kernel we obtained the interacting density 
response function of the g quantity which enabled 
us to calculate the correlation energy via an integration 
over the strength of the Coulomb interaction. The results 
turned out to be rather close to accurate Monte-Carlo re- 
sults. From the response function one can also indirectly 
obtain a local vertex function which can be used to cal- 
culate the electron self-energy and thus the one-electron 
Green function. Unfortunately, this local vertex turned 
out to be too simple minded to remedy the infamous 
band-narrowing proble m^^i^° in simpler metals - and in 
transition metals for that matter.- - In the second paper^^ 
we went beyond the EXX approximation and calculated 
the XC potential at the level of the self-consistent GW 
approximation in atoms. This potential turned out to 
be superior to that of the second-order M0ller-Plesset 
(MP2) approximation,— at least in the outskirts of the 
atoms. The resulting ionization potentials turned out to 
be surprisingly accurate. In the third paper^^ we calcu- 
lated the EXX kernel at imaginary frequencies and used 
it to calculate static polarizabilities, van der-Waals coef- 
ficient and correlation energies of closed-shell atoms. The 
polarizabilities and van der-Waals coefficients were very 
similar to those of time-dependent Hartree-Fock theory 
and thus not very accurate. (See also Ref. H-) On the 
other hand, we obtained better than 95% of the corre- 
lation energies which should be compared to errors of 
the order of a factor of two within the RPA. In contrast, 
our calculated discrete particle-conserving excitation en- 
ergies were not far from those of the RPA. We attributed 
this partial failure to a poor description of the ground 
state within the EXX approximation. 

It is interesting to investigate the performance of the 
TDEXX approximation as far as concerns the descrip- 
tion of photoabsorption spectra. In previous papers we 
were, relying on a new numerical technique based on cu- 
bic splines. While this technique turned out to be ideally 
suited for the necessary inversion of the non-interacting 
density response function and also for the calculation of 
discrete excitation energies, the method was less effec- 
tive in the case of continuous spectra. This deficiency is 
mainly a result of a too sparse and unpredictable sam- 
pling of the continuum especially at larger energies, some- 
thing which would require a prohibitively large number 
of splines to remedy. For these reasons we have in the 
present work designed interpolation schemes and new nu- 
merical methods for directly extracting the important pa- 
rameters determining the shape of the Fano resonances^ 
resulting from the existence of quasi discrete inner-shell 
levels interacting with continua. 

While we have seen that the energy dependence of the 
EXX kernel can improve on many of the deficiencies of 
typical adiabatic (read static) approximations we have in 
the present work made a very important although some- 
what disappointing discovery. The TDEXX completely 
fails to describe parts of the spectra, especially at higher 
energies. At least if the spectra are calculated from 
Fourier transforms to frequency space as most researches 



consider to be the most natural approachi^ Moreover, 
this failure is not connected, e.g., to the neglect of pro- 
cesses describing correlation effects. The failure is caused 
by the existence of zero eigenvalues of the non-interacting 
response function at certain energies. Any method for 
constructing an XC kernel starting from the linearized 
Sham-Schliiter (LSS) equation will suffer from the same 
ailment no matter how much of correlation one is try- 
ing to account for. And those variational approaches to 
many-body theory which build on the Klein functional 
will always result in the LSS equation, suggesting that a 
reasonable theory for constructing better kernels should 
start from a more elaborate variational functional like, 
e.g., that due to Luttinger and Ward. But this is, of 
course, only speculation at the present stage. Such in- 
vestigations will be the topics of future works. 

The paper is organized as follows. In Sec. II we present 
the linear response matrix equation for excitation ener- 
gies and oscillator strengths and discuss how a frequency- 
dependent kernel might influence the results. In Sec. Ill 
we present the EXX equations. We also discuss possible 
generalizations and point out inherent limitations and 
potential problems. In Sec. IV the numerical method is 
briefly discussed. Sec. V is devoted to numerical results 
with a presentation of the EXX kernel and the conse- 
quences of zero eigenvalues in the KS response function. 
The excitation energies and continuum spectra of Be and 
Ne are also presented and compared to the results of the 
RPA and the adiabatic EXX (AEXX) approximations. 
In Sec. VI we give our conclusion and summarize our 
findings and, finally, in the Appendix we derive the Fano- 
profile formula^ using the adiabatic TDDFT formalism. 

II. LINEAR RESPONSE EQUATION 

Within TDDFT the electronic linear density response 
function x is given by 

X ^ Xs + Xs{v + Uc)x, (1) 

where Xs is the Kohn-Sham (KS) linear density response 
function, v is the Coulomb interaction and /xc is the XC 
kernel defined as the functional derivative of the XC po- 
tential Dxc, 



Eq. ([T]) can be considered as a matrix equation in r and 
r' or in any other suitable basis. 

The KS linear density response function or, equiva- 
lently, the retarded KS polarization propagator, can be 
written as 

r,r ,w = > — — — -. 3 

Here, q = (fc,/i) is a composite index, in which k la- 
bels an occupied state and /i an unoccupied state. The 
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excitation function fq = (r)iy9fe(r) is a product of an 
occupied and an unoccupied KS orbital and = — 
is a KS excitation energy. The tilde in fq signifies a mul- 
tiplication by 2yJU^. With Xs in this representation, Eq. 
([T|) for the full x can be rewritten as2^ 

X(r,r',c.) ^ 5]/,(r)[(c. + ^Q-^f\ - Y]-^]f;,{r') (4) 

qq' 



where I is the identity matrix and 

yqq' = ^qSqq' + {fq\v + kc{^)\fq')- 



(5) 



The second term represents an integral over all spatial 
variables of f + /xc multiplied with two excitation func- 
tions fq and /*, . The matrix V is symmetric and if, in 
addition, /xc is frequency-independent and hence real, V 
can be diagonalized as U^VU = 7? , where Z is a diago- 
nal matrix. Then Eq. Q simplifies to 



X(r,r',w) 



(6) 



Here, ^^(r) = Y^q'^qq' fq'{r), where U is the unitary 
matrix which diagonalizes V. The square root of the 
eigenvalues can be interpreted as the new transition fre- 
quencies and cjq as the excitation amplitudeSfSL For most 
known, or generally used, adiabatic kernels the matrix V 
is dominated by the Coulomb interaction and therefore 
its eigenvalues are positive. As a consequence x has the 
correct analytic structure with a positive spectral func- 
tion. It is also easy to see that both x and Xs have the 
same large frequency behavior and thus both obey the 
/-sum rule. 

A frequency-dependent /xc with the correct analytic 
structure requires /xc to have both real and imaginary 
parts. The matrix V is then not Hermitian. It is, 
however, still symmetric and can, in general, be diag- 
onalized according to U'^VU = Z2, where the eigen- 
values and eigenvectors now can be both complex and 
cij-dependent. A particular consequence of having a 
frequency-dependent kernel is that more zeros in the de- 
nominator of Eq. ([6]) can be generated, which would 
imply more resonances in x- These resonances could cor- 
respond to multiple-particle excitations or other excita- 
tions of collective nature. Unfortunately, an approximate 
frequency dependence is not guaranteed to yield a posi- 
tive spectral function. It could also lead to poles in the 
upper half of the complex plane, thus destroying the an- 
alytic structure of x- We will later demonstrate that this 
is indeed a reality to be considered. 



III. EXACT-EXCHANGE EQUATIONS 

The derivation of the EXX potential and kernel has 
been given in several previous publications by ^ ^10, 14^22 
and othersj ^^'^^ Our choice of approach is to start from 



that variational approach to many-body theory which 
builds on the Klein functionalji^ This functional of the 
Green function gives us the total action of the system and 
it contains the functional whose derivative with respect 
to the Green function gives the electronic self energy, 
S = 1^1^ When possible Green functions are restricted 
to those which can be generated by local multiplicative 
potentials the action becomes a functional of the den- 
sity. The stationary property of the action with respect 
to variations in the particle density yields the linearized 
Sham-Schliiter (LSS) equation^ for the KS potential Wxc 
(with ri^i — > 1): 



X.(l, 2)«xc(2)d2 = / S,(2, 3)A(3, 2; I)(i2d3, (7) 



where Ss is the self-energy calculated with KS orbitals 
generated by Wxc and 

zA(3, 2; 1) = ^^^^ = G.(3, 1)G,(1, 2). 

Here, Gg is the KS Green function and V is the total ef- 
fective potential of the KS system. A further variation of 
the LSS equation with respect to the potential V results 
in an equation for /xc: 



X.(l,2)/,,(2,3)x.(3,4)<i2<i3 
^||i2A(3,2;l)<i2.3 

A(1,2;4)A(2,3)G,(3, 1)d2(i3 
Gs(l,2)A(2,3)A(3,l;4)(i2d3, (8) 



where A(2, 3) = S,(2, 3) - Wxc(2)(5(2, 3). 

In the EXX approximation one chooses the HE approx- 
imation for the functional 



$ - -Tr [GGv] , 



(9) 



and the variation of the resulting self energy with respect 
to V becomes 



'^^^(2,3) 
5V{A) 



= -f(2,3)A(2,3;4). 



The terms on the right hand side of Eq. ^ are, in this 
approximation, represented diagrammatically in Fig. [TJ 
The first order vertex diagram, the first on the second 
row in Fig. [1]), which will be referred to as Rv^ is given 
in terms of KS orbitals ipk and eigenvalues Sk as 

Rv{v, r', ^) = -2 ^ ^ (p,, (r)^^^ (r)^*, (r')^fe^ (r') 
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,/x. 





FIG. 1: Diagrammatic representation of the equation for /x. 

where z ~ lu + iO^ and nj. is the occupation number of 
state k. The other four diagrams are self-energy terms 
and their sum, referred to as is: 

kik2k3 

4 [ (rifc, -nfcj(£fc, -£fc3) 



Z2 - (£fc^ - £fcJ2 



(10b) 



We see immediately that the poles of Ry and are 
located at the KS eigenvalue differences, just as in the 
case of Xs- The EXX approximation can, therefore, not 
describe multiple-particle excitations. In order to incor- 
porate those one has to include diagrams which contain 
new poles, as obtained, e.g., by choosing S in the GW 
approximation. The new poles are in that approximation 
generated by the dynamically screened interaction, W . 

It is apparent from Eq. ([1]) that only the quantity 
Xs/xc is needed to obtain the spectrum. Consequently, it 
is sufficient to invert Xs once in Eq. ([5]). It has, however, 
been known for some time that Xs is not always invertible 
due to the existence of nontrivial vanishing eigenvalues 
at particular frequencies >2i At these frequencies there ex- 
ists an external perturbation which does not produce a 
density response to first order. Thus, in the unlikely case 
that R = Rv + -Re has zero eigenvalues at the same fre- 
quencies the kernel /x divergesi^ This can have a drastic 
effect on the calculated spectra as we will see later. 



IV. NUMERICAL APPROACH 

For the numerical implementation we have used an ap- 
proach based on cubic splines as radial basis functions. 
This approach was used also in previous work s "'^^'22 and 
showed to be ideally suited for solving equations where 
an inversion of Xs is needed. In short, a cubic spline 
is a piecewise third order polynomial with compact sup- 
port. Defined on four sub-intervals, and hence composed 
of four different third order polynomials with in total 16 
unknown constants, it is uniquely determined up to a 
multiplicative constant by imposing continuity up to the 
second derivative. A mesh is distributed on the r-axis up 
to a finite Tmax and the basis set is formed by constructing 
a spline starting at every mesh point and extending over 
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FIG. 2: The eigenvalues of Xs and x for Be close the first two 
inner-shell excitations. Two eigenvalues of Xs pass through 
zero. As a consequence, x in TDEXX has zero eigenvalues at 
the same frequencies but with a different oj-behavior. 



four intervals. In this way the splines only overlap with 
its three nearest neighbors on either side. The basis set 
is complete on the interval [0, rmax] in the limit N —^ 00, 
where N is the number of splines. A review on the use of 
general B-splines in electronic structure calculations can 
be found in Ref. ;33'. 

The numerical solution of Eqs. ([7][H]) involves the cal- 
culation of products of KS orbitals. The completeness 
of the spline basis allows us to re-expand these products 
in splines, thus casting the original problem into a linear 
system of equations. The accuracy of this procedure was 
checked by increasing the number of splines and a fast 
convergence was observed. 

For the discrete part of the spectrum in the TDEXX 
approximation we have extracted both excitation ener- 
gies and oscillator strengths from the position and the 
hight of the peaks of Imx- In the RPA and the AEXX 
the kernels are frequency independent and thus excita- 
tion energies and oscillator strengths can be directly ob- 
tained from the eigenvalues and eigenvectors of matrix 
V, defined in Eq. ([5]). We will now describe the method 
that we used in the continuos part of the spectrum. 

Once the radius Tmax is fixed the system studied is ac- 
tually an atom in a box. It turns out, however, that the 
discrete positive energy states provide a good description 
of the true continuum orbitals at the same energy but 
with a different normalization. To get the correct nor- 
malization the box states only need to be multiplied with 
a local density of states factor. By choosing an appro- 
priate mesh both bound and continuum orbitals are well 
described within the same basis. In our calculation we 
have used a mesh which starts as cubic and then smoothly 
changes to a linear mesh, according to the formula 



(11) 
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The photoionization cross sections have been obtained 
by first interpolating the discretized ImXs up to large en- 
ergies, which is possible since it is a smooth function of w 
for every r and r'. The real part of Xs was then obtained 
through the Kramcrs-Kronig relations. This was done for 
every continuum channel. Since the bound states are not 
coupled to the continuum in the KS system, the discrete 
part of Xs can simply be added to form the total non- 
interacting response function of the atom. The resulting 
response function is then used to solve Eq. ^ and to 
obtain the interacting continuum. The accuracy of this 
method has been checked by calculating the RPA spectra 
of some atoms and comparing them with those of other 



V. RESULTS AND DISCUSSION 

In this section we present numerical results on the EXX 
kernel for the Be and Ne atoms. The photoabsorption 
spectra calculated using this kernel as well as its adia- 
batic counterpart are presented and compared. We be- 
gin with a study of the vanishing eigenvalues of the KS 
linear density response function Xs ■ This analysis will be 
useful to interpret and understand the structure of the 
EXX kernel in the frequency domain. 

A. Vanishing eigenvalues of Xs 

Mearns and Kohn,'^^ showed that Xs can have nontriv- 
ial vanishing eigenvalues for finite systems at certain fre- 
quencies. Physically this means that there is a monochro- 
matic perturbation which differs from a constant, and 
which yields a vanishing density response. Note that this 
is true only if the perturbation is not switched on at a 
particular time. If that is the case there is always a finite 
density responsci^ 

The condition for having a vanishing density response 
6n, at frequency loq, can be expressed as 

5n{r) = J2 I dv'f;{v')av') = 0, (12) 

where C is an eigenvector of Xs('j-'o) corresponding to an 
eigenvalue equal to zero and describes the spatial part of 
the monochromatic perturbation. With only one occu- 
pied orbital like, e.g., in He, the /^-functions are linearly 
independent. In that case 5n vanishes only if the inte- 
gral in r' is zero for every g, a condition which implies 
C(r) to be independent of position. On the other hand, 
with more than one occupied level the /^-functions are 
linearly dependent and hence it is, in general, possible to 
find a C, different from a constant, that fulfills Eq. (fT2|) . 

The Be atom has two occupied levels and is there- 
fore the simplest closed-shell spin compensated atom for 
which a nontrivial zero eigenvalue can occur. We have 
diagonalized Xs(w) (in the EXX approximation) for the 



1 = 1 excitation channel of the Be atom in a large finite 
box (40 a.u.) and found a number of eigenvalues which 
pass through zero. The first zero can be found around 

2 a.u., an energy far above the first ionization threshold. 
At higher frequencies, a set of zeros are found close to 
the excitation energies of the inner Is-shell. The eigen- 
values of Xs for frequencies in the range of the Is — > 2p 
and Is 3p transitions are displayed in Fig. [21 We note 
that close to the Is 2p transition (3.948 a.u.) there is 
an eigenvalue which crosses zero at 4.005 a.u., and that 
close to the Is 3p transition (4.058 a.u.) a second 
eigenvalue crosses zero at 4.072 a.u.. We have in fact 
observed that there is one vanishing eigenvalue close to 
every inner-shell transition. 

The Ne atom has three occupied shells and, therefore, 
a larger number of vanishing eigenvalues are expected. 
Indeed, vanishing eigenvalues where found already at low 
frequencies before the first ionization threshold. We did, 
however, not find any vanishing eigenvalues before the 
first excitation energy, in agreement with the findings in 
Ref. [M 

For the numerical calculations of the eigenvalues the 
atom has been confined to a finite box, albeit large. The 
response function is then a sum of discrete transitions and 
it is real almost everywhere. In the limit of an infinite 
box Xs acquires a finite imaginary part at frequencies 
in the continuum. The zero eigenvalues may then also 
acquire a finite imaginary part. 

As seen from Fig. [2] the eigenvalue of Xs which passes 
through zero is proportional to u;^ — ujq for frequencies 
close to loq. The inverse of Xs will, therefore, diverge as a 
simple pole. If the sum Ry + -Rs in Eq. (flO)) has a finite 
component along the eigenvector with zero eigenvalue it 
may have severe consequences for the calculation of /x, 
as we will see shortly. 

A natural way of handling problems associated with 
zero eigenvalues of matrices which need to be inverted 
is so called singular value decomposition. This means 
that one finds the eigenvector corresponding to the zero 
eigenvalue and then invert the matrix in the subspace 
orthogonal to that eigenvector. Notice, however, that 
eigenvalues of Xs are frequency dependent and projecting 
out states within a limited range of frequencies will result 
in abrupt changes in the spectrum at the borders of the 
frequency range. And we can find no physical motivation 
for such drastic manual changes. 



B. The EXX kernel 

As noted previously ; ^^i'^^ the XC kernel is not uniquely 
defined by Eq. ([5]) . Given an f^c one can always add two 
arbitrary functions gi(r,uj) and (72(1"', w) without chang- 
ing the results for x- This originates in the fact that the 
potential is determined up to the addition of a purely 
time-dependent function and that the density variations 
must integrate to zero. When we discuss /x and its de- 
pendence on the frequency we instead consider the quan- 
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FIG. 3: The real an imaginary part of /x for Be around the 
first singularity. 



titles f^"^ (w) = ((7|/x(ti;)|g'), which are unique since the 
excitation functions Integrate to zero and thus remove 
the effect of adding functions of the form gi and 52 • 

Let us first consider the case of the He atom, or for 
that matter, any other spin-compensated two-electron 
system. In this case /x equals minus half the Coulomb 
potential {—■^v). This can be seen either by noting that 
the EXX potential becomes —^fnv and hence the func- 
tional derivative can be taken explicitly, or after some 
manipulations of the diagrammatic expression in Eq. 
([To)) . In the latter case It is Important to Impose the 
LSS equation as well as the additional constrainij^^ 

{^PnIvx — SJ|<pAr) =0 if and only if 11m Vx{r) = 0, 

r — >QO 

where N denotes the highest occupied orbital. We thus 
see that for He the EXX kernel is energy Independent. 

Let us now consider systems with more than one oc- 
cupied orbital. The Be atom has two filled closed shells 
- the inner Is-shell and the outer 2s-shell - and a full nu- 
merical solution is therefore required. The kernel for the 
I = 1 channel was calculated with a finite box radius of 40 
a.u.. At low frequencies /x has a weak frequency depen- 
dence and is purely real. At larger frequencies, however, 
/x develops a strong frequency dependence and around 2 
a.u. the kernel diverges. In Fig. [3] we show the real and 
imaginary parts of /x around this point. The divergence 
is typical of a double pole with the real part given by 
the difference between a term 1/(lu — wq)^ and a squared 
delta function, S{uj — i^-'o)^- A small asymmetry around 
the pole can be observed which Indicates the presence of 
an additional simple pole. The exact location of the pole 
can be traced back to the first vanishing eigenvalue of Xs 
(see discussion in the previous section). The pole struc- 
ture of /x shows that R — Ry + Rt. does not compensate 
for the zero of Xs ■ 

Continuing to higher frequencies more singularities in 
the kernel are observed. Near the Is 2p transition 
in Xs a second pole In /x develops. As a matter of fact 
we find a pole in /x close to every Is np transition 



in Xs (see the upper panel of Fig. 4(a) I. These poles 
have the same structure as the one at 2 a.u. and are also 
confirmed to originate from the vanishing eigenvalues of 
Xs- The vanishing eigenvalues corresponding to the first 
two inner-shell excitations are displayed in Fig. [2l The 
double and simple pole structure in /x is expected since, 
in order to calculate the kernel, Xs has to be Inverted 
twice, see Eq. (|5]). Let us call C,i{r,Lu) the eigenvector of 
Xs(r, r', uj) with eigenvalue Then from the equation 

XsfxXs = R the kernel can be expressed as 



Ciir,Lu)C{r',u;). (13) 



If we denote by the frequency for which there is an 
eigenvalue 1^ equal to zero, we have that Ik ~ — uj'^ for 
UJ ^ ujk and Eq. psp can be cast in a more transparent 
form 



k [ 



(^2 



,2^2 



(14) 



where the quantities , fl and fl are weakly depen- 
dent on w, and the r and r' dependence has been sup- 
pressed. Notice that the object Xs/x, the basic quantity 
needed to calculate the full response function x, only has 
simple poles at the same frequencies. 

The kernel for the Ne atom has an even larger number 
of poles due to the larger number of inner-shells and thus 
a larger number of vanishing eigenvalues of its Xs- AH 
poles are of the same structure as those discussed for 
Be, but the first pole appears already between the first 
and second resonance In Xs: i-e., below the ionization 
threshold of Ne. 

As can be seen from Eq. the /^-representation 

is not well suited for studying the singular behavior of 
the kernel In real space. The elements of the matrix 
V in Eq. (O will all have the pole structure described 
above. When the matrix is diagonalizcd, however, this 
pole structure survives in only one of the eigenvalues, 
whereas the remaining ones appear to have a smooth fre- 
quency dependence but still affected by the pole giving 
rise to an oscillatory behavior. In Figs. 4(a) and |5(b)| two 
different smooth eigenvalues of this matrix are shown for 
Be and Ne respectively. The eigenvalue corresponding 
to the 2s —^ 2p transition in Be is weakly dependent on 
the frequency and hence the intersection with the line lo 
giving the location of the new excitation energy is only 
slightly shifted from that of the adiabatic approximation. 
The new position in the TDEXX approximation is in bet- 
ter agreement with the corresponding excitation energy 
of the TDHF approximation as compared to the same po- 
sition in the adiabatic case, see Tab. IV Bl The position 
of the first discrete excitation in Ne is, however, strongly 
modified by the frequency dependence as is shown in Fig. 
5(a) The kernel of the adiabatic approximation moves 



the peak to a lower frequency as compared to the RPA 
while the TDHF result lies at a higher frequency (0.674 
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FIG. 4: (a) Two different eigenvalues of the matrix V for Be. The upper/lower panel shows an eigenvalue corresponding to 
the first inner/outer-shell excitation. The intersection with the dashed line gives the positions of the new resonances. For the 
inner-shell eigenvalue this occurs exactly where /x has a pole. At that value the oscillator strength is, however, zero. Two more 
solutions are found but for complex uj. (b) The imaginary part of x is plotted on the complex plane. Close to every inner-shell 
transition two poles are found in the complex plane. 




FIG. 5: (a) The absorption spectrum of Ne around the 2p 3s transition. The adiabatic approximation to /x moves the RPA 
value in the wrong direction since the real part of /x is negative at w = 0. The strong frequency dependence of /x moves the 
peak closer to the TDHF values as it should. Also a small negative peak is generated close to the pole of /x at 0.73 a.u.. (b) 
The lower panel shows the eigenvalue corresponding to the 2p — > 3s transition. Three intersections are found. Two of these 
are seen in the left figure while the third has a zero oscillator strength. 



a.u.). The frequency dependence of the kernel leads to 
the opposite effect and moves the peak even a little be- 
yond the true TDHF value. At a somewhat larger fre- 
quency the EXX kernel has a pole and the eigenvalue 
corresponding to the 2p —> 3s excitation is seen to os- 
cillate (see lower panel of Fig. |5(b)[ ). The maximum of 
the oscillation occurs exactly at the location of the pole. 
Apart from the desired solution, i.e., the crossing with 
the Lu-line at 0.680 a.u., there are two more crossings. 



The first of these corresponds to the small negative peak 
found in the spectrum (at 0.71 a.u.) and the next has 
zero oscillator strength. Clearly, the negative peak is un- 
physical but we expect that the ability of the TDEXX to 
describe the spectrum rapidly decreases as we approach 
the pole. 

Since there is a pole close to every inner-shell excitation 
the eigenvalues oscillate in these regions. The eigenvalue 
corresponding to an inner-shell excitation gives only one 



8 



solution exactly at the pole as seen in the upper panels 
of Figs. 



4(a) and 5(b) 



This is true for all inner-shell 
transitions in Be and Ne and seems to be a general be- 
havior. The oscillator strength is zero at these points 
leading to a complete disappearance of the peaks from 
the spectrum as also seen, e.g., at third crossing corre- 
sponding to the first excitation in Ne as discussed above. 
Because the sum rule is obeyed the missing oscillator 
strength is expected to be transferred to excitations of 
the outer shell. This conjecture is based on having the 
correct analytic structure of the response function with 
a Lehmann representation leading to the fact that the 
large frequency behavior of the response is determined 
by the sum of the different oscillator strengths. Unfor- 
tunately, in many cases the three solutions at the real 
frequency axis close to the double pole structure, which 
we discussed above, instead leads to one solution at the 
real axis and two more complex solutions one of which is 
in the upper half plane and the other in the lower. Thus 
the proper analytic structure of the response is destroyed 
and the fact that the TDEXX response has the correct 
large frequency behavior does not guarantee that the full 
spectrum contains the correct sum of oscillator strengths. 
Part of the spectrum is actually missing which is a clear 
failure of the theory. In Fig. |4(b)| x clearly exhibits 
two peaks located symmetrically on opposite sides of the 
real axis. Apart from the solution with zero oscillator 
strength there are thus two more solutions at complex 
frequencies {u). 

Finally, before further discussing the spectrum we will 
make a small digression and investigate the eigenvalues 
of X- In Fig- m the eigenvalues of both Xs and x are 
presented in a frequency range where there are two van- 
ishing eigenvalues of Xs- We see, however, that also x 
have vanishing eigenvalues at the same points. We also 
see that the eigenvalues of x only touches the zero in- 
tensity axis whereas the eigenvalues of Xs cross this axis. 
This can be explained by studying 



X = 



Xs 



1 



VXs 



Xs'R 



(15) 



obtained by using the equation /x = xJ^^Xj^- Our cal- 
culations show that the matrix R does not have a zero 
eigenvalue where Xs has. Consequently, by writing this 
equation in the basis of the eigenvectors of Xs, we im- 
mediately see that x behaves as (w^ — ujq)'^ close to ujq. 
There is no reason to believe that the Coulomb interac- 
tion would change the simple zeros of Xs to double zeros. 
We therefore expect also x to have simple zeros although 
shifted relative to those of Xs ■ From the definition of the 
exact XC kernel. 



/x 



X 



Xs 



V, 



(16) 



/xc should have simple poles where Xs and x have vanish- 
ing eigenvalues. The double pole structure found within 
the EXX approximation is thus an artifact of this theory. 



C. Discrete excitation energies 

Before the first ionization threshold the atomic pho- 
toabsorption spectra consists of a set of discrete tran- 
sitions. In Tab. IV Bl the first two discrete transitions 
are presented for Be and Ne. A comparison is made be- 
tween the AEXX, TDEXX, TDHF and the RPA as well 
as the KS eigenvalue differences of the exact^ and the 
EXX potentials. The latter potential was also used when 
calculating x in the aforementioned approximations. 

The TDEXX approximation is expected to give results 
close to TDHF. In this sense the outer-shell transitions 
of Be are well described already in the AEXX approx- 
imation. By also accounting for the frequency depen- 
dence the results improve even further. For the outer- 
shell excitations of Ne dynamic effects in /x arc crucial. 
The AEXX approximation reduces the values of the RPA 
whereas the TDHF results are larger than those of the 
RPA. The frequency dependence corrects this tendency 
yielding results in good agreement with TDHF, although 
somewhat over-estimated. The effect of the kernel on 
the first excitation energy of Ne is illustrated in Fig. 
5(a)[ The results of TDHF and hence of TDEXX dif- 
fer markedly from the experimental values. By looking 
at the exact KS eigenvalue differences it can be concluded 
that the error is mainly related to the effect of the ground 
state potential. The kernel has only a small effect on the 
outer shell excitation energies as was previously pointed 
out by Petersilka et. al..— 

A description of inner-shell excitation energies in terms 
of density functional eigenvalues is bound to fail. In 
Tab. IVBI we see that the true core excitation energies are 
largely under-estimated by even the exact KS eigenvalue 
differences as compared to the true excitation energies. 
An adiabatic kernel is unable to correct these errors as 
confirmed by the examples (RPA and AEXX) shown in 
the table. A strong frequency dependence in the kernel 
might, however, improve the situation considerably. The 
vanishing eigenvalues of Xs result in a strong frequency 
dependence in the EXX kernel near these excitations. 
As we saw in Sec. IIIC this frequency dependence is, 
however, too strong and even produces the wrong ana- 
lytic structure of the response. The matrix V contains 
separate diagonal blocks for different excitation channels 
and off diagonal blocks for the coupling between them. 
Notice, however, that this separation is not entirely well 
defined since part of the inter-channel coupling is already 
contained in /xc. By including /x and diagonalizing only 
the block containing the inner-shell excitations we obtain 
some improvement for the position of these. Accordingly, 
for Be we have diagonalized the Is np block and for Ne 
we have diagonalized the 2s — > np block. The results are 
marked c in Tab. IV Bl The most important observation 
here is the fact that by means of this somewhat ad-hoc 
procedure we can obtain results for those inner-shell ex- 
citations which completely vanish in the full treatment. 

The calculation of the inner-shell transitions was done 
with a finite box size and the excitation energies where 
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Trans. 


KS(EXX) 


RPA 


AEXX 


TDEXX 


TDHF" 


KS (Exact)' 


Exp.'' 


ls-^2p 
ls^3p 
2s^2p 
2s^3p 


3.948 
4.058 
0.131 
0.241 


3.959 
4.060 
0.203 
0.255 


3.956 
4.059 
0.177 
0.247 


4.005" 
4.072" 
0.176 
0.247 


4.346 
4.646 
0.176 
0.247 


4.017 
4.153 
0.133 
0.269 


4.243 
4.461 
0.194 
0.274 


Ne 

2s^3p 
2s^4p 
2p^3s 
2p^4s 


1.604 
1.667 
0.659 
0.779 


1.608 
1.668 
0.667 
0.781 


1.607 
1.667 
0.663 
0.780 


1.609" 
1.669" 
0.680 
0.783 


0.674 
0.782 


1.542 
1.602 
0.612 
0.725 


1.674 
1.731 
0.619 
0.727 



"From Ref. IliandliO 
''From Ref. |4l| 

'^Transition calculated assuming no intor-shell coupling. 
''From Ref. El 



TABLE L The two first discrete excitation energies from the 2s and Is shell of Be and the 2p and 2s shell of Ne. Different 
approximations for the kernel is used in conjunction with the EXX potential for the ground state. A column with the eigenvalue 
differences calculated using the exact KS potential of Umrigar et. aL*^ is also presented. All values are in a.u.. 
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FIG. 6: (a) Photoionization cross section for Be after the first ionization threshold, 
from the Is 2p and Is 3p transitions. 



(b) First two Fano resonances resulting 



extracted from the location of the delta peaks in Imx. In 
reality these peaks will be somewhat shifted and broad- 
ened forming a so-called Fano resonance due to the cou- 
pling to continuum channels. In the next section we will 
present results for these resonance structures. 



D. Photoionization cross sections and 
Fano resonances 

The 2s — > 2p photo ionization cross section of Be is 
presented in Fig. 6(a) and different approximations are 
compared. The experimental result has not been in- 
cluded since in this region the spectrum is dominated by 
the 2p, ns and 2p, nd double excitation resonances, which 
cannot be captured by the approximations studied here. 



The figure shows that the AEXX and the TDEXX ap- 
proximations are almost indistinguishable and also very 
close to the RPAE results of Amusia et. al.— We remind 
the reader that the RPAE approximation is identical to 
linearized TDHF. The frequency dependence of the ker- 
nel is weak up to around 2 a.u. where the first singularity 
of /x occurs. This strong frequency dependence of /x has, 
however, only a small effect on the spectrum in the stud- 
ied region. 

The photoionization cross section after the first ioniza- 
tion threshold of Ne is presented in Fig 
two ionization channels, 2p 



7(a) 



It contains 
to continuum s or d. The 
AEXX agrees very well with the RPAE results.-^ There 
is also a fairly good agreement with experiment. In the 
figure we have deliberately not presented results obtained 
by including the full energy dependence of /x. In the 
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energy region covered by Fig. 7(a) the non-interacting 



response function of Ne has a number of zero eigenval- 
ues producing violent structures in f^- This unphysical 
behavior results in a series of strong peaks also in the 
full response. Thus including the energy dependence of 
/x results in absorption spectra with no resemblance to 
the spectra in Fig. 7(a)[ These strong structures are of 
course obtained using discrete splines which give only a 
discrete sampling of the continuum. It could be possible 
that calculating /x using true continuum functions would 
somehow alleviate this problem even though it will not 
go away completely. The zero eigenvalue below the ion- 
ization threshold does produce unphysical structure in 
the absorption spectra of Ne. In this case one cannot 
blame the unphysical spectrum on the lack of a proper 
continuum, and we are convinced of the futility in go- 
ing through the trouble of calculating proper continuum 
functions. 

At higher frequencies, Fano resonances due to excita- 
tions of inner-shell electrons occur. In the Appendix we 
present a derivation of the Fano profile formula for excita- 
tions with single-particle character starting from the lin- 
ear density response function within adiabatic TDDFT. 
The Fano parameters g, F, and are there given ex- 
pressions in terms of the adiabatic kernel /xc, see Eq. 
(|A7l) and Eq. ([ATT|) . In Figs. [6(b)] and [7(b)] we compare 
spectra calculated either directly from Eq. ^ using an 
interpolated Xs or by employing the analytic method de- 
scribed in the Appendix. Within the analytic approach 
we have for both Be and Ne ignored the coupling be- 
tween different discrete states via the continuum (which 
corresponds to ignoring the off-diagonal elements of the 
matrix M in Eq. (|A12P ). For Ne we have also assumed 
that we can ignore the coupling between the s and d con- 
tinua, which can be reached from the 2p-shell. The Figs. 
[6(b)[ and [7(b)[ contain results of both the RPA and the 
AEXX approximation. As seen, the two different ways 
of calculating the spectra around a resonance give almost 
identical results. This verifies the validity of the assump- 
tions made above. 

The two first inner-shell resonances of Be are very 
sharp with a zero on one side of the resonance due to 
the coupling to only one continuum. The very large 
g-parameter is accompanied with a very small width, 
q sa 220 and F w 1.6 meV in the RPA and q w 320 and 
rwO.86 meV in the AEXX approximation, for the case of 
the Is 2p resonance. In Ref. [471 the same parameters 
were calculated using the ALDA kernel giving similar re- 
sults for the g-parameter. The width (F) from the AEXX 
approximation is, however, a factor of two smaller than 
those from both the RPA and the ALDA. Experimental 
results have in this case not been found. 

In Ne the discrete excitations of the 2s-shell couples to 
both the p and the d continua. The minima on one side of 
the resonances therefore never reaches zero, which means 
that < 1. This parameter is for the first resonance 
p2 w 0.17 in the RPA and ~ o.60 in the AEXX approx- 
imation. Since the experimental value is 0.70 the AEXX 



approximation improves over the RPA. The widths are 
F « 15.4 meV in the RPA and F « 14.4 meV in the 
AEXX approximation. Again, the value of the AEXX 
approximation is in quite good agreement with the ex- 
perimental result of 13 meV. The q parameter differs, 
however, substantially from the experimental values. The 
RPA gives 6.3 and the AEXX approxi mation 3.5 whereas 
the experimental value is 1.6. In Fig. [7(b)[ the two first 
resonances of Ne are displayed. We clearly see that the 
shape is very sensitive to the g-parameter. Notice that 
the experimental^ and the RPAE results of Amusia^^ are 
shifted in order to facilitate the comparisons. The RPAE 
approximation yields q — 0.95 and F = 34 meV. The 
latter is about three times larger than the experimental 
value. In conclusion wc can say that for the 2p — > 3s reso- 
nance of Ne neither of the static approximations here give 
a reasonable account of the experimental findings. The 
deviations to experiment is in many cases much larger 
than the results from different theories. A more sophis- 
ticated theory is therefore required. Theories which do 
provide a reasonable description of the spectra are wave 
function based, include many terms in a configuration 
interaction expansion, and can hardly be used on larger 
systems. 



VI. SUMMARY AND CONCLUSIONS 

In the present work we have used the so called exact- 
exchange approximation within TDDFT to calculate par- 
ticle conserving discrete excitation energies and photoab- 
sorption spectra of a few spherical atoms. As in sev- 
eral earlier papers on this topic we have also here re- 
lied on a numerical method based on cubic splines as 
basis functions for one-electron wave functions, polariz- 
abilities, and the XC kernel. Clearly, such a numerical 
method amounts to a discretization of the continuum in 
which we have no control over the resulting positions of 
one-electron eigenvalues or electron-hole excitation en- 
ergies which are ingredients in the necessary correlation 
functions. In order to circumvent such problems we have 
here constructed interpolation schemes which work very 
well for the description of, e.g., the non-interacting KS 
response function and other functions in the continuous 
part of the spectrum. In addition, we present in the Ap- 
pendix a new way in which the parameters of the Fano 
profiles associated with the autoionizing resonances cor- 
responding to inner-shell excitations can be calculated 
directly within adiabatic TDDFT without resorting to 
a full calculation of the entire spectrum. The resulting 
parameterized spectra are actually very accurate approx- 
imations to the full spectra which was demonstrated for 
the RPA and for the spatially non-local EXX kernel eval- 
uated at zero frequency (AEXX). But, unfortunately, it 
has not been true when we have tried to incorporate the 
full frequency dependence of the EXX kernel. For rea- 
sons which have been discussed in the paper and will 
be mentioned again later on, the violence of that fre- 
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FIG. 7: (a) Photoionization cross section for Ne just after the first threshold, (b) First two Fano resonances arising from the 
2s — > 3p and 2s — > 4p transitions. Notice that the experimental and the RPAE results have been shifted to allow for a better 
comparison of widths and g-parameters. 



quency dependence has defied any reasonable interpola- 
tion scheme. As a result we present few results based on 
the fully frequency dependent EXX kernel. 

One such result is the calculation of discrete transi- 
tion energies. There are really no strong trends in the 
results. We here compare our results to the correspond- 
ing results from a full time-dependent Hartree-Fock (HP) 
calculation. As discussed in the paper, this is motivated 
by our belief that the latter calculation is the 'target' for 
the TDEXX approximation. It is not motivated by the 
HF results being very close to experiment which they are 
not. We find that the rather poor eigenvalue differences 
in the EXX approximation are improved by the inclu- 
sion of the Hartree part of the kernel (the RPA results). 
Then, for Ne, the situation deteriorates when the static 
(adiabatic) exchange effects are included just to again im- 
prove when the exchange effect are treated in a frequency 
dependent fashion. The full frequency-dependent kernel 
has, however, a tendency to overcorrect the errors from 
the RPA. In the case of Be the results of the adiabatic 
approximation instead represent a marked improvement 
on the RPA results and the frequency dependence gives 
a slight additional improvement. In all fairness we must 
add that these feeble trends are completely overshadowed 
by starting from a better XC potential like the exact DF 
potential available for the light atoms or from the very 
similar XC potential from a full GW calculation. The 
outer transition energies then move much closer to ex- 
perimental results. 

In the continuous part of the spectra we have used our 
new way of obtaining the parameters of the Fano line 
shapes and compared our results to those of the RPA, 
to previous calculations using different versions of the 
ALDA, and to experiment. We have found that all adi- 
abatic kernels including the AEXX give a clear improve- 
ment over the RPA results for those Fano parameters 



which, roughly, correspond to the weights and the widths 
of the resonances with a slight edge for the AEXX ap- 
proximation. The fact that the parameter mainly respon- 
sible for the height of the resonances is poorly described 
by all static approximations suggests that the former re- 
sults are fortuitous. This impression is further reenforced 
by observing that a supposedly superior treatment in- 
cluding full energy dependence like the RPAE yields re- 
sults for the parameters which are further away from ex- 
periment. Turning then to the density- functional version 
of RPAE, i.e., the TDEXX the description completely 
breaks down and the resonances disappear as will now 
be discussed. 

Perhaps the most important finding of the present 
work is the violent frequency variation of the EXX ker- 
nel discovered in the vicinity of energies associated with 
inner-shell excitations. This, in turn, is a result of van- 
ishing eigenvalues of the non-interacting Kohn-Sham re- 
sponse function close to these energies, and the necessity 
to invert that response function in order to obtain the re- 
sulting EXX kernel - at least in the way we presently for- 
mulate the theory. Although we can observe a tendency 
of the EXX kernel to use a strong frequency variation as a 
way to mimic the effect of the non-local but frequency in- 
dependent three-point vertex of time-dependent Hartree- 
Fock theory this frequency behavior is quite unphysical 
and actually leads to the complete disappearance from 
the spectrum of all inner-shell excitations. We here trace 
this breakdown of the TDEXX approximation to the ap- 
pearance of poles at a finite distance in frequency into the 
upper complex frequency plane thus producing a density 
response function with an incorrect analytic structure. 
One might think that this could be an order of limits 
issue and that the problem might go away if a proper 
continuum was included in the calculations. We doubt 
the correctness of this conjecture because the problem is 
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already visible in the truly discrete part of, e.g., the Ne 
spectrum. In the present approach we implicitly calculate 
the linear response to a perturbation which has a sinu- 
soidal variation at all times. It has been shown by van 
Leeuwen^^ that the non-interacting response function is 
always invertible for perturbations which completely van- 
ish before a given time. This suggests that a proper 
TDEXX spectrum could be obtained by applying a si- 
nusoidal perturbation at one particular time, propagate 
the time-dependent EXX orbitals, wait a time T until all 
transients have died out and then Fourier transform the 
charge density starting from time T to a very large time. 
This would amount to some limiting procedure which is 
quite feasible to carry out. It is, however, hard to see 
that such an approach would lead to something different 
to what we already have obtained. 

It should be noted that the inclusion of higher-order 
correlation effects like, e.g., screening or relaxation ef- 
fects will not be a remedy to the problem. Any proce- 
dure for obtaining better XC kernels based on the vari- 
ational formulation of many-body theory and the Klein 
functional will lead to the linearized Sham-Schliiter equa- 
tion at the first variation, albeit with a very sophisticated 
self-energy containing many physically important corre- 
lation effects. One further variation of the LSS equation 
in order to obtain the kernel will again reveal the neces- 
sity to invert the non-interacting Kohn-Sham response 
and the problem is still there. Only by pure coincidence 
will the right hand side of the equation determining the 
kernel have a zero eigenvalue at the same frequency as 
the non-interacting response. A possible resolution to 
the problem could be starting from a more sophisticated 
and variationally stable functional than the Klein func- 
tional, e.g^ the functional of Luttingcr and Ward (LW), 
see Ref. |32|. The additional computational work will be 
substantial already at the exchange-only level but it is 
definitely worth a serious research effort. But then, of 
course, we are no longer working within the TDEXX. 
The underlying ^-functional will still be at the level of 
the HE approximation but it could be hoped that the 
more complicated LW functional would render a kernel 
within TDDET which would result in a response function 
much closer to that of TDHE. The latter is expected to 
have the correct analytic structure but it does not pro- 
duce overly impressive results. 
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APPENDIX A: FANO RESONANCES 

In this section we derive expressions for the Fano 
parameters^ of the autoionizing resonances of single- 
particle character. The derivation is based on adiabatic 
linear response TDDFT. 

In the zeroth-order approximation, given here by the 
non-interacting KS system, there is no coupling between 
single-particle transitions. Different excitation channels 
can therefore be treated independently. When interac- 
tions are present and two channels overlap, one with dis- 
crete levels and the other with continuum levels, the cou- 
pling forms resonance structures having so-called Fano 
profiles. 

Starting from the assumption of having only one dis- 
crete transition superimposed on only one continuum, the 
KS response function can be written as 

Choosing the retarded response function which is ana- 
lytic in the upper half of the complex plane we have 
z = u + iO~^ . The quantities /o and Eq £^re the ex- 
citation function and the excitation energy of the dis- 
crete transition respectively. Notice that we have defined 
/o(r) = 2y^/o in order to make the formulas as light as 
possible. The ionization potential / is naturally smaller 
than £o- As an example of this model we can consider 
the Is — > 2p transition of Be where the other Is — > np 
transitions are neglected and the continuum is given by 
2s ep. 

Within TDDET the fully interacting can be ob- 
tained from Eq. (|3|): 

X^(r,r',z) = /o(r)[z2l_V]oo7o(r') 

+JdsUr)[zH~YU'fo{r') 

+jde' Uv)[zH-Y]-^,fAv') 

+ jdede' Uv)[zH-\]-Jh,{v'). (A2) 

The matrix V has a submatrix containing elements in- 
volving only excitations to the continuum. This matrix 
can be diagonalized initially by a linear transformation of 
all the continuum excitation functions — J de'U^^i f^i . 
If the kernel is frequency independent (adiabatic) this 
transformation is unitary. The matrix z^I — V then takes 
the form 

f Vof, Voe' . . . \ 

Veo {z^ - e^)S,,, ^ (A3) 

where Vqo = £q = £q + voo, and T^o = Vqs = -voe- 
From Sec. IH vqo = {fo\v + /^d/o) and vq, = {fo\v + 
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/^l^e), where the superscript A signifies an adiabatic 
approximation to /xc- 

Analytic expressions for the elements of the inverse of 
z^I— V can now easily be obtained. Defining the complex 
function 

i?(r,z)= fds^, (A4) 



J — £2 



we can write a-s 



X'^(r,r',z) = 



(/o(r) + S(r,z))(/o(r') + i?(r',z)) 



el-F{z) 



(A5) 



_ r°°ric 9e (r)ge (r ) ^ j ^j^g non-rcsonant back- 



where Xs — Ji de 
ground continuum and 



(A6) 



From now on, let us assume that we have integrated 
X^ with some perturbing time-dependent potential and 
thereby removed the dependence on r and r'. After some 
manipulations we can extract the real and imaginary 
parts of x^ ■ If w > we have 



ReF = V I de . °" ImF = -7:-°^ 



3 
2uj 



Note that the real and imaginary parts of F depend on 
oj. For w > we also have that 



2uj 



If we define 



-2Imi^, e 



fo + RcB 
Im B 



{lu^ -el-ReF) 



(A7) 



we can write 



+ 1 

The full x^ can then be written as 



Imx^ + Rex^ (A8a) 
(A8b) 



X = 

[e + i) 



Imx^ + Rex^ (A9a) 



If < the first term in Eq. (jA8a[l acquires a minus 
sign and the full x^ becomes 

X^(-) ^^^^^^Imx^ + Rex^ (A9b) 



These equations show that x^ has the correct analytic 
structure. It is analytic in the upper half of the complex 
plane and has a symmetric real part and an antisymmet- 
ric imaginary part. If we redefine F — > F/{lu + Eq) such 
that ReF ReF/{uj + eo) and F r/(tj + eo) the poles 
of x^ are located at So + Re FztiT /2, where the plus sign 
refers to < 0. With a reasonably well behaved kernel 
the poles are in the lower half of the complex plane as 
they should. The real part corresponds to the position 
of the resonance and the imaginary part to the width. 

The formula in Eq. ljA8b[) gives the asymmetric res- 
onance profile due to Fano but the parameters are here 
defined in terms of TDDFT quantities. 

With this derivation we have shown that Fano res- 
onances for single-particle excitations occur already in 
RPA, i.e., with /xc = 0, and how an adiabatic kernel 
modifies the parameters. (See also Ref. |49| ) 

All parameters are frequency dependent, but this is 
assumed to be weak and approximately constant over 
the resonance region. Double and multiple-particle ex- 
citations are not present in the adiabatic approximation 
since such states are absent in the KS response function. 
Maitra et. ali^ have studied how /xc should behave in 
order to describe double excitations. 

The above discussion has been limited to the case of 
one discrete state and one continuum. The formula can, 
however, easily be generalized to the case of several dis- 
crete states coupled to several continua. The derivation 
does not involve steps more complicated than those above 
and therefore we simply state the formulas and discuss 
the parameters. If we have n discrete states and m con- 
tinua the response function can be generalized to {uj > 0) 



Rex«(c.) = 



^'^^;.f^^% ^Imxg + Rexg 

i'^t+Qr? 2 



Imx» = 5:^^J^p?Imx? + (l-^p,2)I^^) 

where Xs = Sl'Li Xs fc' i-^-j a sum over the non-resonant 
background for every continuum channel. We see that 
we have different g-parameters for every discrete state 
as well as different widths. These are given by the new 
definitions of F and B which now become Fi — ^2^=1^!^ 
and Bi = X^fcLi ^f- The new parameter p is defined as 

r,/2 Imx^ 

and gives the fraction of the continuum which mixes with 
the discrete state. If there is only one continuum p = 1 
for every discrete state and we have a zero on one side of 
the resonance. In order to obtain Eq. (|A10|1 for x''^ we 
have assumed that the matrix 



(A12) 



where 



k=i 
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is diagonal. Physically, the matrix W describes the inter- 
action between the discrete states via the continua and 
the neglect of ofF-diagonal elements of this matrix has a 
very small effect on the spectra as seen by comparing our 
parametrized spectra to the full ones in Sec. IV Dl 



The derivation above is based on adiabatic TDDFT. If 
the kernel is frequency dependent it also has an imaginary 
part which cannot be neglected. We notice that with a 
general complex kernel it is not possible to derive the 
equations above using the same steps. 
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